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£S) ' Abstract. 

5_i ■ Is gravitational growth responsible for the observed large scale struc- 

ture in the universe? Do we need non-gaussian initial conditions or non- 
gravitational physics to explain the large scale features traced by galaxy 
surveys? I will briefly revise the basic ideas of non-linear perturbation 
theory (PT) as a tool to understand structure formation, in particular 
through the study of higher order statistics, like the skewness and the 
3-point function. Contrary to what happens with the second order statis- 
tics (the variance or power-spectrum), this test of gravitational instability 
■ is independent of the overall amplitude of fluctuations and of cosmic evo- 

lution, so that it does not require comparing the clustering at different 
CO ' redshifts. Predictions from weakly non-linear PT have been compared 

•O with observations to place constraints on our assumptions about struc- 

ture formation, the initial conditions and how galaxies trace the mass. 
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1. Introduction 



Where does structure in the Universe come from? The current paradigm is 
that it comes from gravitational growth of some small initial fluctuations. The 
self-gravity of an initial overdensed region increases its density contrast so that 
eventually the region collapses. For a flat Universe in the linear regime, the 
local density contrast 5 = p/p — 1 grows as the expansion factor, eg D = a, so 
that since decoupling linear gravitational growth has the potential of amplifying 
fluctuations by at least a factor of a thousand. But Gravity is not linear and 
when objects start collapsing the growth could be much larger. On galactic 
scales one also has to consider other forces such as hydrodynamics, heating and 
cooling by friction, dissipation, feedback mechanism from stars, such as nova 
and supernova explosions, interaction with the CMB and so on. 

To test if the above picture of gravitational growth is correct we need to 
deal with a classical initial condition problem. Because gravitational time scales 
are very slow, we have no way to measure the growth of individual large scale 
structures and we need to resort to the statistical study of mean quantities. One 
can imagine, for example, measuring the rms fluctuations (at a given scale) at 
different cosmic times to see if this agrees with the predicted amount of grav- 
itational growth, D. Observationally this corresponds to finding the clustering 
properties of some tracer of structure (eg galaxies) at different redshifts. If the 
tracer is not perfect, we will have some statistical biasing. The problem with 
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this approach is that by the time the rms fluctuations change significantly there 
typically has also been a substantial cosmic evolution of the corresponding trac- 
ers. Thus, it is difficult to disentangle the effects of the underlying cosmological 
model (which sets the rate D of gravitational growth) from galaxy evolution. 

It is therefore important to have a way of testing the gravitational growth 
paradigm at a single cosmic time or redshift. Higher order correlations and 
weakly non-linear clustering allows us to do just this. This is because one can 
construct ratios of higher order correlations to powers of the two point ampli- 
tude which are independent of cosmic time or cosmological parameters, but still 
contain information of the underlying dynamics. 

2. Gravitational Growth in the weakly non-linear regime 

Gravitational growth increases the density contrast of initially small fluctuations 
so that eventually the region collapses. The details of this collapse depends on 
the initial density profile. As an illustration we will focus in the spherically 
symmetric case. Thus we will study structure growth in the context of matter 
domination, the fluid limit and the shear-free or spherical collapse approxima- 
tion. These turns out to be very good approximation for the one-point cumulants 
of the density fluctuations. It is then easy to find (see Peebles 1980, Gaztanaga 
& Lobo 2000) the following second order differential equation for the density 
contrast 6 in the Einstein-deSitter universe (Qk = Q\ = 0), eg a(t) = (t/to) 2 ^ 3 : 

d 2 5 1 dS 3 4 1 (d5\ 2 3 2 

d 2 V + 2 dq ~ 2 3 l + <5 \dq) + 2 U 

where we have shifted to the rhs all non-linear terms, and used the rj = ln(a) 
as our time variable. This equation reproduces the equation of the spherical 
collapse model (SC). As one would expect, this yields a local evolution so that 
the evolved field at a point is just given by a (non-linear) transformation of 
the initial field at the same point, with independence of the surroundings. The 
linear solution factorizes the spatial and temporal part: 

Si{x, t) = D(t - t ) S(x,to) = D(t) 5 (x) (2) 
where D is the linear growth factor, which from the above differential equation: 

D(t) = Ci e" + C 2 e- 3 / 2r ? = d a(t) + C 2 a(t)- 3 / 2 (3) 

with growing D ~ a and the decaying D ~ a~ 3 / 2 modes. Thus, the initial 
fluctuations, no matter of what amplitude, grow by the same factor, D, and the 
statistical properties of the initial field are just linearly scaled. For example the 
linear rms fluctuations 07 or its variance a 2 gives: 

a 2 = £ 2 = (5 2 (t)) = (D(t-t ) 2 6 2 ) = D(t-to) 2 a 2 (4) 

where Sq = 6 (to) and ao refer to some initial reference time to: a 2 =< 5q >. 

We are interested in the perturbative regime (5 — > 0), which is the relevant 
one for the description of structure formation on large scales. The non-linear 
solution for 5 can then be expressed directly in terms of the linear one, 5f. 
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n=l 

Thus all non-linear information is encoded in the v n coefficients. We can now 
introduce this expansion in our non-linear differential equation, with Si given by 
the linear growth factor D = a = e v , and compare order by order to find: 

34 682 446440 8546480 

V2 = 21 ; ^ 3 = l89 ; VA = 13659" ; ^ = M3243" (6) 

These results are derived for the Einstein-de Sitter, but are also a good approxi- 
mation for other cosmologies (eg Bouchet et al. 1992, Bernardeau 1994a, Fosalba 
k, Gaztahaga 1998b, Kamionkowski & Buchalter 1999). For non-standard cos- 
mologies or a different equation of state see Gaztahaga &; Lobo (2000). 

One can now find the N-order cumulants £ N = (S N } c , where N = 2 cor- 
responds to the variance. Here the expectation values (...) correspond to an 
average over realizations of the initial field. On comparing with observations 
we assume the fair sample hypothesis (§30 Peebles 1980), by which we can com- 
mute spatial integrals with expectation values. Thus, in practice (...) is the 
average over positions in the survey area. It is useful to introduce the N-order 

hierarchical coefficients: SV = > e § skewness for 53 and kurtosis for £4. 

These can easily be estimated from the series expansion above by just taking 
expectation values of different powers of 8 (eg see Fosalba & Gaztahaga 1998a). 
For leading order Gaussian initial conditions we have: 

S 3 = 3v 2 ; S 4 = 4 u 3 + 12 v\ ; 5 5 = 5 v 4 + 60 u 3 u 2 + 60 v\ (7) 

These results have also been extended to the non-Gaussian case, see Fry h 
Scherrer (1994) Chodorowski &; Bouchet (1996), Gaztahaga & Mahonen (1996), 
Gaztahaga & Fosalba (1998). If we take for v 2 the non-linear solution above, 
eg v 2 = 34/21, the skewness yields S3 = 3^2 = 34/7, which reproduces the 
exact perturbation theory (PT) result by Peebles (1980). Thus the above (SC) 
model gives the exact leading order result for the skewness. This is also true for 
higher orders (see Bernardeau 1992 and Fosalba Sz Gaztahaga 1998a). These 
expressions have to be corrected for smoothing effects (see Juszkiewicz et al. 
1993, Bernardeau 1994a, 1994b, and Fosalba & Gaztahaga 1998a) and possibly 
from redshift space distortions (eg Hivon et al 1995, Scoccimarro, Couchman and 
Frieman 1999). Next to leading order terms have been estimated by Scoccimarro 
& Frieman (1996), Fosalba & Gaztahaga (1998a, b). 

The 1-point cumulants measured in galaxy catalogues have been compared 
with these PT predictions (eg Bouchet et al. 1993, Gaztahaga 1992, Gaztahaga 
& Frieman 1994, Baugh, Gaztahaga & Efstathiou 1995, Baugh & Gaztahaga 
1996, Colombi etal 1997, Hui & Gaztahaga 1999). The left panel of Figure 1 
shows a comparison of Sn measure in the APM Galaxy Survey (Maddox et al. 
1990), with the predictions above (see Gaztahaga 1994, 1995 for more details). 
The agreement between predictions (lines) and measurements (points) on scales 
R > 10/i -1 Mpc (where £ 2 < 1)) is quite good indicating that the APM galaxies 
follow the non-linear gravitational growth picture. For errors on statistics see 
Szapudi, Colombi and Bernardeau (1999) and references therein. 
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Figure 1. Figures in the left panel show estimates of Sjy, for = 
3 — 6, from the APM galaxy catalogue. Lines correspond to the PT 
predictions. Figures in the right panel shows the APM projected 3- 
point amplitude 53(0), for triplets of objects 1, 2, 3 on the sky at fixed 
#12 = #13 = 2 deg and as a function of a the interior angle between 
these two triangle sides. The thick line shows the PT predictions. 



2.1. Biasing: tracing the mass 

The expressions above apply to unbiased tracers of the density field; since galax- 
ies of different morphologies are known to have different clustering properties, 
at least some galaxy species are biased. As an example, suppose the probability 
of forming a luminous galaxy depends only on the underlying mean density field 
in its immediate vicinity. The relation between the density field as traced by 
galaxies S ga i(x) and the mass density field 5(x), can then be written as: 

«*) = /(<*(*)) = £ ^ S n ( x ), (8) 

n U - 

where b n are the bias parameters. Thus, note how biasing and gravity could 
produce comparable non- linear effects. To leading order in £ 2 , this local bias 

scheme implies £2" = ^1^2 an< ^ ( see Fry & Gaztanaga 1993): 

Q gal S3 o ^ 2 Q90.1 ^4 . -, ^gjj , h n &2 / n \ 

5 3 =-^ + 3^ , 5 4 =^ + 12-^ + 4^ + 12^ , ... (9) 

Gaztanaga &: Frieman (1994) have used the comparison of S3 and S4 in PT with 
the corresponding measured APM values (as shown in Figure 1) to infer that 
b\ ~ 1, 62 — and 63 ~ 0, but the results are degenerate due to the relative 
scale-independence of Sn and the increasing number of biasing parameters. One 
could break this degeneracy by using the configuration-dependence of the pro- 
jected 3-point function, 93(a), as proposed by Frieman & Gaztanaga (1994), Fry 
(1994), Matarrese, Verde & Heavens (1997) Scoccimarro et al (1998). As shown 
in Frieman & Gaztanaga (1999), the configuration-dependence of 93(a) on large 
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scales in the APM catalog is quite close to that expected in perturbation theory , 
suggesting again that b\ is of order unity (and 62 — 0) for these galaxies. This is 
illustrated in the right panel of Figure 1. The solid curves show the predictions 
of weakly non-linear gravitational growth. The APM galaxy measurements are 
shown as symbols ; other curves show results for each of the zones. The agree- 
ment indicates that large-scale structure is driven by non-linear gravitational 
instability and that APM galaxies are relatively unbiased tracers of the mass on 
these large scales. 

3. Conclusions 

The values of Sn = 5,^/5,2 1 can be measured as traced by the large scale 
galaxy distribution (eg Bouchet et al. 1993, Gaztahaga 1992, 1994, Szapudi 
el at 1995, Hui & Gaztahaga 1999 and references therein), and also the weak- 
lensing (Bernardeau et al. 1997, Gaztahaga & Bernardeau 1998) or the Ly- 
alpha QSO absorptions (Gaztahaga & Croft 1999). These measurements of the 
skewness S3, kurtosis £4, and so on, can be compared with the predictions from 
weakly non- linear perturbation theory (see Figure 1) to place constraints on 
our assumptions about gravitational growth, initial conditions or biasing at a 
given redshift (see Mo, Jing & White 1997). Contrary to what happens with the 
second order statistics (eg the variance), this test of gravitational instability is 
quite independent of the overall amplitude of fluctuations and other assumptions 
of our model for cosmological evolution, and does not require comparing the 
clustering at different redshifts. As shown in Gaztahaga & Lobo (2000), one can 
also use the Sjv measurements to constraint non-standard cosmologies. 

Frieman & Gaztahaga (1999) have presented new results for the angular 
3-point galaxy correlation function in the APM Galaxy Survey and its compar- 
ison with theoretical expectations (see also Fry 1984, Scoccimarro et al. 1998, 
Buchalter, Jaffe & Kamionkowski 2000). For the first time, these measurements 
extend to sufficiently large scales to probe the weakly non-linear regime (see 
previous work by Groth & Peebles 1977, Fry & Peebles 1978, Fry & Seldner 
1982). On large scales, the results are in good agreement with the predictions 
of non-linear perturbation theory, for a model with initially Gaussian fluctua- 
tions (see Figure 1). This reinforce the conclusion that large-scale structure is 
driven by non-linear gravitational instability and that APM galaxies are rela- 
tively unbiased tracers of the mass on large scales; they also provide stringent 
constraints upon models with non-Gaussian initial conditions (eg see Gaztahaga 
& Mahonen 1996; Peebles 1999a,b; White 1999; Scoccimarro 2000). 
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